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Abstract. In this paper we characterize the slack matrices of cones and poly- 
topes among all nonnegative matrices. This leads to an algorithm for deciding 
whether a given matrix is a slack matrix. The underlying decision problem is 
equivalent to the polyhedral verification problem whose complexity is unknown. 



1. Introduction 

This paper is concerned with a class of nonnegative matrices with real entries, 
called slack matrices, that arise naturally from polyhedral cones and polytopes. 
Given a polytope P C K" with vertices V\, . . . ,v p and facet inequalities ajx < f3j 
for j = 1, . . . , q, a slack matrix of P is the p x q nonnegative matrix whose (i,j)- 
cntry is f3j — ajvi, the slack (distance from equality), of the ith vertex Vi in the jth 
facet inequality ajx < /3j of P. A similar definition holds for polyhedral cones. 

Slack matrices form an interesting class of nonnegative matrices with many spe- 
cial properties. Most obviously, if M is a slack matrix of a polytope P, then the 
zeros in M record the face lattice of P and hence the combinatorial structure of 
P. In its entirety, M specifies an embedding of P up to affine transformation. 
However, slack matrices carry much more surprising information about P. In [15j . 
Yannakakis proved that the nonnegative rank of a slack matrix of P is the minimum 
k such that P is the linear image of an affine slice of the positive orthant K/L. We 
use R+ to denote the set of nonnegative real numbers. The nonnegative rank of a 
matrix M £ R^* 9 is the smallest k such there there exists vectors 01, . . . , a p £ 
and bi, . . . , bg £ R*L such that Mij — afbj. Affine slices of positive orthants that 
project onto P are called polyhedral lifts or polyhedral extended formulations of P 
and the smallest k such that admits a lift of P is called the (polyhedral) ex- 
tension complexity or nonnegative rank of P. If the extension complexity of P is 
small (polynomial in the dimension of P), then usually it is possible to optimize a 
linear function over P in polynomial time by optimizing an appropriate function on 
the lift. This is a powerful technique in optimization that yields polynomial time 
algorithms for linear optimization over complicated polytopes. There are many 
instances of n-dimensional polytopes with exponentially many (in n) facets that 
allow small polyhedral lifts. 
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Yannakakis' result was generalized in [5] to lifts of convex sets by affine slices of 
convex cones via cone factorizations of slack operators. Even in the larger context 
of cone lifts of convex sets, the case of polytopes is the simplest and the theory 
relies on slack matrices of polytopes and their factorizations through cones. Thus, 
understanding the structure of these matrices is fundamental for this theory. There 
are several phenomena that occur in the class of nonnegative matrices that have 
not yet been observed for slack matrices. For instance, an important open question 
is whether there exists a family of slack matrices of polytopes that exhibit an 
exponential gap between nonnegative rank and positive semidefinite rank. (If Si_ 
denotes the cone of k x k real symmetric positive semidefinite matrices, then the 
positive semidefinite rank of a matrix M £ R^* 9 is the smallest k such that there 
exists matrices Ai £ Si, i = 1, . . . ,p and Bj £ S^_, j = 1, . . . , q such that = 
(Ai, Bj}.) While there are simple families of matrices that exhibit even arbitrarily 
large gaps between nonnegative and positive semidefinite ranks [5] Example 5], no 
family of slack matrices with this property is known. Such a family would be a 
clear witness for the power of semidefinite programming over linear programming 
in lifts of polytopes. 

This paper was motivated by the many open questions about slack matrices 
which rely on understanding the structure of these matrices. We establish two 
main characterizations of slack matrices of polyhedral cones and polytopes. In 
Section [2] we establish linear algebraic characterizations: Theorem [1] for cones and 
Theorem |6] for polytopes. In Section |4] we give combinatorial characterizations: 
Theorem |2"21 for polytopes and Theorem [2^1 for polyhedral cones. In Section [3] we 
use our characterization from Section [2] to give an algorithm for recognizing slack 
matrices. The computational complexity of this problem is unknown and is equiva- 
lent to the polyhedral verification problem. There are several further geometric and 
complexity results about slack matrices throughout the paper. 

Notation: For a set of vectors A — {ai, . . . , a p }, cone(^l) := XiOi : \ > 0} 
is the cone spanned by A; conv(^l) := {^2 A.^a.; : Ai > 0, — 1} is the convex 
hull of A; lin(_4) := \ai : A; € M} is the linear span of A, and aff(_4) := 
Xidi : A; = 1} is the affine span of A. The above sets can also be defined for 
an infinite subset A C 1™ by taking unions over all finite subsets of A. For a n x q 
matrix M, we let rows(M) and cols(M) denote the sets of all rows and columns, 
respectively, of M. We let A ■ M be the set of vectors {x T M : x £ .4}. For a set 
K C R™, lineal(i^) is the largest subspace contained in K, known as the lineality 
space of K . The dimension of a polytope P, dim(P) is the dimension of aff(P), the 
affine hull of P, and the dimension of a cone K is the dimension of Ym(K). 

2. Geometric Characterizations of Slack Matrices 

2.1. Slack Matrices of Polyhedral Cones. Consider the polyhedral cone 

K = {x £ R n : x T B >®} = W + -A 

in R n constrained by the columns of the matrix B £ W lXq and generated by the 
rows of the matrix A £ R px ™. We call (the set of rows of) A a V -representation and 
(the set of columns of) B an H-representation of K. The slack matrix of K with 
respect to the representation (A, B) is S = AB £ R^ 9 . Its (i, j)-entry records the 
"slack" of the ith generator of K with respect to the jth inequality of K in the 
given description of K. 
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Let Sk denote the set of all slack matrices of K. For S £ Sk, any matrix 
obtained by scaling the rows and columns of S by positive reals is again in Sk 
since scaling the vectors in a V and/or ^-representation of K does not change 
K . Also, Sk can have matrices of different sizes as adding redundant inequalities 
and/or generators to the representations of K does not change K. From 

(R" • B) n R q + = K ■ B = (R p ■ A) ■ B = R p ■ S 

C (R p • s) n R% = (R p • AB) ni« C (R" • B) n R+ 

we find that R p + • S = R p ■ S PI R+ which says that the cone generated by the rows 
of S coincides with the nonnegative part of the row span of S. In fact, this relation 
characterizes slack matrices of cones: 

Theorem 1. A nonnegative matrix M £ R pX9 is a slack matrix of a polyhedral 
cone if and only if 

(1) R p + ■ M = R P ■ M DR q + , 

or in other words, the cone spanned by the rows of M coincides with the nonnegative 
part of the row span of M . 

Proof. It remains to show that every matrix M £ R p + Kq with R^ • M = W ■ M C\M.\ 
is a slack matrix of some cone. Let n = rank(M) and choose a bijective linear map 

tp : MP ■ M -4 R" 

that preserves the (standard) scalar product (an isometry). Let M, denote the ith 
row of M and let A £ R px " be the matrix whose rows are ip(Mi). Let ir : R 9 — >• 
R p ■ M be an orthogonal projection and let B £ W ixq be the matrix whose columns 
are y>(7r(<Bi)), . . . (p(ir(<B q )) where <Ej is the ith standard unit vector in R g . Then 
M = AB and using ©, 

K = {x£ R" : x T B >€)} = {ip(y) : y G R p • M, tp{y) J B > O} 

= ip(W ■ M n R^_) = yj(R+ • M) = R p • A , 
which shows that M is a slack matrix of the cone K. □ 
Recall that the dual cone of K is the cone 

K* = {y £ R n : x T y > for all sc G K] = {y £ R™ : Ay > 0} 
Hence, S 111 is a slack matrix of if* and we get the following result. 

Proposition 2. A nonnegative real matrix is a slack matrix of a polyhedral cone 
if and only if its transpose is also the slack matrix of a polyhedral cone. 

In particular, we obtain the following consequence of Theorem [1] 

Corollary 3. A nonnegative matrix M £ R px<3 is a slack matrix of a polyhedral 
cone if and only if 

(2) M-R q + =M-R«nR p , 

or in other words, the cone spanned by the columns of M coincides with the non- 
negative part of the column span of M. 

We say that a matrix M satisfies the row cone generating condition (RCGC) 
if (p} holds and the column cone generating condition {CCGC) if J2]) holds. 
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Corollary 4. For a nonnegative matrix M G l^j* 9 the following statements are 
pairwise equivalent: 

• M is a slack matrix of a polyhedral cone. 

• M satisfies the RCGC. 

• M satisfies the CCGC. 

The equivalence of RCGC and CCGC for a general nonnegative matrix is not 
obvious. However, its proof becomes transparent via the theory of slack matrices 
of polyhedral cones and cone duality. 

For a nonnegative M with RCGC/CCGC, the proof of Theorem Q] showed how 
to produce a cone K such that M G Sk ■ We now give another way to produce such 
a cone K that will be useful later. For any matrix M G R pxq of rank k, we will call 
a factorization of the form M = AB with A G M? xk , rank(A) = k and B G R kxq , 
rank(-B) = k a rank factorization of M. 

Lemma 5. Let M G R^ X9 be the slack matrix of a polyhedral cone and let M = AB 
be a rank factorization of M . Then if K is the cone generated by the rows of A, 
the columns of B form an Ti-representation of K . In particular, M G Sk- 

Proof. Let K = cone({ai, . . . , a p }) and K = {x E R k : x T V > 0, j = 1, . . . , q} 
where a, is the ith row of A and b> is the jth column of B. We need to show that 
K = K. Since M is a nonnegative matrix, K C K. To prove the reverse inclusion 
we will argue that every linear function that is nonnegative on K is also nonnegative 
on K. Let k = rank(M). Since M has the CCGC and A ■ R fc = M ■ R q (since 
rank(A) = fc), we have that M-M.% = A-R k r)R+. Suppose L(x) = l\X\+. . .+4x fc > 
for all x G K. Then the evaluation vector (L(ai), . . . , L(a p )) T — Ai lies in 
M ■ R\. Since the columns of M are AW for j = 1, . . . , q, there exists Xj > such 
that At = \j(AV) = ^Ej=i x j bj ■ This implies that I = Y, q j=1 XjV since 

the columns of A are linearly independent. Therefore, i is a nonnegative linear 
combination of the 6 J 's and L(x) > is valid on K . □ 

2.2. Slack Matrices of Polytopes. We now investigate the slack matrices of 
polytopes. Let V G R pxn and P = conv(rows(F)) be the polytope in W 1 that is 
the convex hull of the rows of V. Suppose also that P = {x G M™ : Wx < w} 
with W G R qxn and w G R q . To avoid unnecessary inconveniences, we assume 
that dim(P) > 1. We call (the set of rows of) V a V- representation and (the set 
of columns of) [w, — W] T an ^-representation of P. The slack matrix of P with 
respect to the representation (V, W, w) is then 

(3) S = [1, V] • [w, —W] T G R p i xq . 

We denote the set of all slack matrices of P by Sp. Clearly, scaling the columns 
of a slack matrix of P by positive scalars yields another slack matrix of P, because 
scaling the vectors in an ^-representation of P yields another H-rcprcscntation of 
P. However, we cannot scale the rows of a matrix S G Sp and still stay in Sp. 

The matrix S is also the slack matrix of the homogenization cone of P: 

(4) P h = R\ ■ [1, V] = {(x Q ,x) elxf : Wx < x Q w} 

r t i 

with respect to the representation ™ w t ). Since dim(P) > 1, there is 

some c G M." with 

max{c T a; : x G P} — min{c T x : x G P} = 1 , 
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and hence, due to LP-duality, we get 

(5) (1, O t ) G R 9 ■ (w, W) and so also, (1, O t ) G R 9 ■ (w, -W). 

From d3j) and (© we get that 1 G S ■ R 9 , the column span of S. These properties 
characterize the slack matrices of polytopes of dimension at least one: 

Theorem 6. A matrix M G R^ X9 with rank(M) > 2 is a slack matrix of a polytope 
if and only if M is a slack matrix of a polyhedral cone and 1 G M ■ R 9 . 

Proof. It suffices to show that a matrix M G R^ 9 with 1 e M that is the slack 
matrix of some cone K C R™ with respect to a representation (A, B) is also the 
slack matrix of some polytope. To construct such a polytope, choose any /iel' 
such that 1 = Mfi and define c = B[i. Then Ac = 1 since M — AB. Define 
P = conv(rows(yl)). Then we have: 

P = {y T A : y T l = l,yeR p + } = {y T A : y T Ac = 1, y G R p + } 

= {x G K : x T c = 1} = {x G R" : x T B > O, x T c = 1} . 

Mapping the hyperplane in R™ defined by x T c — 1 isometrically to R n_1 (as in 
the proof of Theorem [T|) , we find that M is a slack matrix of the resulting image 
of P. □ 

Corollary 7. A matrix M G R+ x<3 with rank(M) > 2 is a slack matrix of some 
polytope if and only if it satisfies the RCGC (or, equivalently, the CCGC) and 
1 E M ■ R 9 holds. 



Theorem [T] geometrically characterizes the slack matrices of cones as those ma- 
trices M G M^ X9 that satisfy 

(6) cone(rows(Af)) = lin(rows(M)) n M.% . 

There is an analogous geometric characterization of slack matrices of polytopes. 

Corollary 8. A matrix M G with rank(M) > 2 is a slack matrix of some 

polytope if and only if 

(7) conv(rows(A/)) = aff(rows(A/)) n R q + . 

Proof. First, suppose that M is a slack matrix of some polytope. Then by Corol- 
lary [3 we have that M satisfies ([6]) and 1 G M -R 9 . Hence, there exists some c G R 9 
such that Mc = 1 and the affine hyperplane L — {x G R 9 : x T c — 1} contains the 
rows of M. Intersecting L with both sides of (O, we obtain ([7]). 

For the reverse implication, let M G R^* 9 be a nonnegative matrix satisfying (J7]). 
Using any isometry ip between the d-dimensional affine subspace aff(rows(M)) 
and R d , we find that M is a slack matrix of the ^-image of the polytope defined 
in 0. □ 

We have seen above that every slack matrix of a polytope P has the all-ones 
vector in its column span and is also a slack matrix of the homogenization cone P h 
of P. The next example shows that not all slack matrices of P h are slack matrices 
of P, in fact, this does not even hold for the slack matrices of P h that have the 
all-ones vector in their column span. 

5 



J. Gouveia, R. Grappe, V. Kaibel, K. Pashkovich, R.Z. Robinson & R.R. Thomas 



Example 9. Let P be the square [— 1, l] 2 . The matrix 



M 



I 



\ 



f 1 


2 2 

3 3 


\ 


1 


1 -1 




1 


-1 1 




V 2 


-2 -2 






is in Sph and 1 is in the column span of M. It is clear, however, that M is not in 
Sp since each facet of [—1, l] 2 is equidistant from the two vertices not on the facet. 
On the other hand, since M has the RCGC/CCGC and 1 is in its column span, 
it is the slack matrix of some other polytope Q. To obtain it, write a new rank 
factorization of M (note that rank(M) = 3) so that the first factor contains the all 
ones vector as its first column as follows: 
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-2 / 
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-2 

















[ 1/4 


1 







V 1/4 
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1 1 


1 
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3/4 -5/4 


-1/4 




1/4 


-1/4 -1/4 


3/4 




5/4 



M 



to get 



M = 



By Lemma [5j M is the slack matrix of the cone with V-representation the rows of 
the first factor and ^-representation the columns of the second factor. Assuming 
the coordinates of this three-dimensional cone are xq,x\,X2, and slicing the cone 
with the hyperplane {((xq, x\, x-i) : xq — 1} gives a polytope Q with vertices 
(2/3, 2/3), (1, —1), (—1, 1), (—2, —2) and ^-representation given by the columns of 
the second factor. Then M € Sq. 

2.3. Further Results on Slack Matrices of Cones and Polytopes. In this 
section we derive some more insight into the geometric relations between cones, 
polytopes, and their slack matrices that will be useful in later parts of the paper. Wc 
return to the setup used earlier: K is assumed to be a cone and S the slack matrix 
of K with respect to its representation (A, B) where A e W xn and B e 1™ X? . 

First, we will show that every slack matrix of a cone is the slack matrix of 
some pointed cone. Recall that we use Ym(K) to denote the linear hull of K and 
lineal^) to denote the lineality space of K. Then we have lin(_ftT) = W ■ A and 
lineal(X) = leftkernel(S). A cone K is pointed if lineal^) = {O}. Define 

L := lm(K) n lineal^)- 1 = (W p ■ A) n (B ■ R q ) . 

Then we have 

\in(K) = L + lineal(K) 

(where the summands are orthogonal to each other) and 

K = (K n L) + lineal(if) , 

where ifHLCIisa pointed (i.e., having trivial lineality space) cone with dim(A"n 
L) = dim(L). Denoting by A' e M. pxn the matrix obtained from A by orthogonal 
projections of all rows to L, we have 



K f) L 



■A' 



and 



A'B. 
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By mapping L isometrically to R dlm ( L ), we thus find that S is a slack matrix of the 
pointed cone that is the image of K l~l L under that map and we get the following: 

Lemma 10. A matrix is a slack matrix of a polyhedral cone if and only if it is a 
slack matrix of some pointed polyhedral cone. 

If the cone K is pointed, then for every zero-row of S — AB the corresponding 
row of A is a zero-row as well. Hence, removing any zero-row from S results in 
another slack matrix of K . A similar statement clearly holds for adding zero-rows. 

Lemma 11. If a matrix S is a slack matrix of a pointed polyhedral cone K then 
every matrix obtained from S by adding or removing zero-rows is a slack matrix 
of K as well. 

Lemmas [TU1 and [TTI together also imply this statement: 

Lemma 12. If a matrix is a slack matrix of some polyhedral cone then every 
matrix obtained from it by adding or removing zero-rows is a slack matrix of some 
polyhedral cone as well. 

Let us further investigate the linear map x i— > x T B. It induces the isomorphism 

(8) L ► W ■ S 

isomorphism 

between the linear space L and the row span of S because of the relations: 
L C lineal^)- 1 = leftkerne^S)^ 

and 

L-B = (L + lineal(iT)) ■ B = ]m(K) ■ B = (M p ■ A) ■ B = W P ■ S . 
It also induces the isomorphism 

K n L ► R p + ■ S 

isomorphism 

between the cone K D L and the cone spanned by the rows of S since (K n L) ■ B = 
{{K n L) + lineal(if)) • B = K ■ B = (R p + ■ A) ■ B = R p + ■ S. In particular, we have 
shown the following result: 

Lemma 13. A polyhedral cone K is pointed if and only if dim(K) — rank(S') for 
any slack matrix S of K . 

Recall that if P is a polytope with representation (V, W, w) and slack matrix 
S= [1,V] B where 



B = 



T 
W 



-W T 

then the homogenization P h of P is a pointed cone that also has S as a slack 
matrix. Since P h is pointed, L contains the entire cone and we can restrict the 
isomorphism in © to the set {1} x P = conv(rows([l, V})). Thus we have that 
{1} x P is isomorphic to conv(rows([l, V])) ■ B — conv(rows(S*)). This establishes 
the first part of the following: 

Theorem 14. If S is a slack matrix of the polytope P, then P is isomorphic to 
conv(rows(S')). In addition, we have dim(P) = rank(S') — 1. 

Proof. To prove the second statement, note that dim(P' 1 ) = dim(P) + 1. By 
Lemma H3l we have that dim(P' 1 ) = rank(iS). □ 
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In the conic case, we had that M £ Sk if and only if M T £ Sk* ■ This correspon- 
dence breaks down for polytopes as we see in the example below. The reason behind 
this is that we cannot scale V-representations of polytopes by positive scalars. 

Example 15. The matrix 



M = 



( 1 


i 








\ 


1 





1 








1 








1 








1 








1 








1 





1 


V o 








1 


1 / 



is a slack matrix for the triangular prism in R 3 . Thus, by Corollary [71 M satisfies 
both the RCGC and the CCGC, and the all ones- vector is in the column span of M. 
However, the all-ones vector is not in the row span of M, so M T is not the slack 
matrix of any polytope. 

Despite this complication, we can still derive some results for transposes of slack 
matrices of polytopes. Recall that the polar of a polytope P C M. n is 

P° = {y£ E™ : x T y < 1 for all x £ P} . 

Then P° is a polytope whenever £ int(P), the interior of P. Since translating 
P does not change its slack matrices, we may assume that £ int(P). Therefore, 
P has an ^-representation of the form P = {x G M™ : Wx < 1} and P° = 
conv(rows(W0). Similarly, if P = conv(rows(V r )), then P° = {x £ R™ : Vx < 1}. 
This implies that the slack matrix of P with respect to the representation (V, W, 1) 
is the transpose of the slack matrix of P° with respect to the representation (W, V, 1) 
and we get the following result that is analogous to Proposition [5] for cones. 

Proposition 16. For any polytope P, there exists a slack matrix M £ Sp such 
that M T is also a slack matrix of a polytope. 

In the light of Theorem El this says that slack matrices of polytopes (which 
already have 1 in their column span) allow positive scalings of their columns that 
puts 1 into their row span as well. This is false for general nonnegative matrices. 

Example 17. Continuing Example [T5l we see that the following matrix M' ob- 
tained by scaling the columns of M is also a slack matrix of the same prism and 
does have 1 in its row span: 



M' = 



( 2 


2 








\ 


2 





4 








2 








4 








2 








2 








4 





2 










4 


2 J 



The prism has vertices: 

(0, 1, -1), (2, -1, -1), (-2, -1, -1), (0, 1, 1), (2, -1, 1), (-2, -1, 1) 

and M' comes from the facet description: 

z < 1, — y < 1,—x + y < 1, x + y < 1, — z < 1. 
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Therefore, P° has vertices (0, 0, 1), (0, -1, 0), (-1, 1, 0), (1, 1, 0), (0, 0, -1) and is a 
bisimplex with slack matrix M lT . 

We can also show a converse to Proposition \W[ 

Proposition 18. Suppose M £ R^ 9 such that M and M T are both slack matrices 
of poly topes. Then there exists a polytope P, with £ int(P), such that M £ Sp 
and M T £ Spo . 

Proof. Since M T is a slack matrix of a polytope, we have that 1 £ M?_ -M. Without 
loss of generality, we can scale M by a positive scalar so that 1 £ conv(rows(M)). 

Let M be a slack matrix of a polytope R with dim(i?) = d. By Theorem IT4l 
rank(M) = d + 1. Since the convex hull of the rows of M is isomorphic to R, 
we have that the affine hull of the rows of M has dimension d. Let J denote the 
all-ones matrix of dimension p x q. Since 1 is contained in the affine hull of the 
rows of M, we have that the affine hull of the rows of M — J passes through the 
origin and has dimension d. Hence, rank(M — J) = d. This implies that we can 
write M-J = AB with A £ W xd and B £ R dxq . 

Let A' = (1, A) and let B' = (l, B T ) T . Then M = A'B' is a rank factorization 
of M. Let P := conv(rows(A)) and Q := {x £ R d : 1 + x T B > ©}. Then the rows 
of A' form a V-representation of P h and the columns of B' form a ^-representation 
for Q h = {(x ,x) e R d+1 : lx + x T B > O}. By Lemma El P h = Q h which 
implies that P = Q. Therefore, M is a slack matrix of P and M T is a slack matrix 
off . " □ 

3. An Algorithm to Recognize Slack Matrices 

In this section, we discuss the algorithmic problem of deciding whether a given 
nonnegative matrix has the RCGC (or, equivalently, the CCGC). According to 
Corollaries |4] and [7] this is the crucial step to be performed in order to decide 
whether a given matrix is a slack matrix of a cone or a polytope. 

We start with a promising result: 

Theorem 19. The problem to decide whether a nonnegative matrix satisfies the 
RCGC (or the CCGC) is in coNP. In particular, the same holds for checking the 
property of being a slack matrix (of a cone or of a polytope). 

Proof. If the given matrix M £ E^ xg does not satisfy the RCGC, then there is some 
point x £ M. p ■ M n M. q + \ K5. • M (which can be chosen to have coordinates whose 
encoding lengths are bounded polynomially in the encoding length of M). The fact 
that x $ K.5. ■ M can be certified by the help of some separating hyperplane whose 
normal vector can be chosen to have coordinates with encoding length bounded 
polynomially in the encoding length of M as well. □ 

Next, we are going to describe an algorithm to check the CCGC (equivalently, 
the RCGC) for a nonnegative matrix. By Corollary 21 this algorithm will then 
provide a method to check if a given nonnegative matrix is a slack matrix of a cone. 
To check if the matrix is the slack matrix of a polytope (see Corollary [7]), we can 
add the additional step of checking if the all-ones vector is in the column span of 
the matrix which is doable in polynomial time. A SAGE worksheet implementing 
this code can be found at http://www.matri.washington.edu/~rzr. 

Algorithm to check if a nonnegative matrix has the CCGC 
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Input: A matrix M £ R^* 9 . 

Output: True if M has the CCGC and False otherwise. 

(1) Compute a basis L for the left kernel of M. For each vector £ in L, generate 
the equation £ T x — 0. 

(2) Generate an ^-representation of the cone K with the equations from the 
previous step and the inequalities X\ > 0, . . . , x p > 0. 

(3) Compute a minimal V-representation of K. 

(4) Normalize the vectors in the V-representation and the columns of M. 

(5) Check that each normalized vector in the V-representation is a normalized 
column of M. If so, return True. If not, return False. 

Proof. We have K = M ■ R q n R p + and M ■ R\ C K due to the nonnegativity of M. 
Thus, M satisfies the CCGC if and only if K C M ■ M. q + holds, which is what the 
algorithm checks in the last three steps (note that all cones involved are pointed 
because they are contained in K+). □ 

The only computationally challenging part of the algorithm is converting from 
the ^-representation of K to a V-representation. There are several algorithms to 
do this, and we refer to jB], [5], and [TT] for information on the different techniques. 
No polynomial time algorithm for this conversion exists, since the V-representation 
may have size exponential in that of the ^-representation. If the dimension of 
the cone is fixed, however, then there do exist polynomial time algorithms for the 
conversion [3]. Thus, we obtain the following complexity results. 

Theorem 20. For fixed r, checking whether a rank r matrix satisfies the RCGC 
(CCGC) can be done in polynomial time. In particular, checking whether matrices 
of fixed rank are slack matrices of cones or polytopes can be done in polynomial 
time. 

Given an %-polyhedron P and a V-polytope Q contained in P, the problem 
of deciding whether P = Q is known as the polyhedral verification problem. The 
complexity of this problem is unknown |13j . However, a polynomial time algorithm 
for the polyhedral verification problem would yield an output sensitive algorithm 
for the problem of computing the facets of a polytope given in V-representation, 
and thus solve a decades old open problem in computational geometry (see [7]). 

Clearly, given a V-polytope it is easy to check whether it is contained in an 
"H-polyhedron. The reverse problem of checking whether an H-polyhedron is con- 
tained in a V-polytope is known to be coNP-complete [4] . Note that the polyhedral 
verification problem is the restriction of the latter problem to those instances in 
which the V-polytope is contained in the %-polyhedron (see also 
http://www.inf.ethz.ch/personal/fukudak/polyfaq/node21.html, [8 and (T2"]). 

Theorem 21. The following problems can be reduced in polynomial time to each 
other: 

(1) The polyhedral verification problem 

(2) Is a given matrix a slack matrix of a polytope? 

(3) Is a given matrix a slack matrix of a cone? 

(4) Does a given matrix satisfy the RCGC /CCGC? 

Proof. Corollary [7] shows that (2) can be reduced (in polynomial time) to (4) (since 
checking whether 1 is contained in the column space can be done in polynomial 
time) and Corollary [4] shows that (4) can be reduced to (3) . 
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We can also reduce (3) to (2): Suppose we need to check whether a given ma- 
trix M is a slack matrix of a cone. By Lemma [JTJ we can assume that M has no 
zero rows. We can also scale the rows of M by positive scalars without effect on M 
being a slack matrix of a cone. Using these two facts, we can assume that 1 is in 
the column span of M. Then, being a slack matrix of a cone is equivalent to being 
a slack matrix of a polytope due to Theorem [SJ 

Since Corollary [5] shows how to reduce (2) to (1), it thus remains to establish a 
reduction of (1) to (2). Let Q = conv(rows(V)) with V G W xn and P = {x G K n : 
Wx < w} with W G R qxn and w e W with Q C P. Suppose we need to decide 
whether P = Q. First, we check whether P is pointed (i.e., W has a trivial right 
kernel) and dim(P) = dim(Q) (both checks can be done in polynomial time, the 
second one using linear programming). If either check fails, then P ^ Q. 

So let us assume dim(P) = dim(Q) and that P is pointed. The latter fact implies 
that the affine map p : M. n — >• R 9 defined via ip(x) = w — Wx is injective. Let M 
be the matrix arising from V by applying ip to each row. Then, due to Q C P, 
we have that M is nonnegative. According to Corollary the matrix M is a slack 
matrix of a polytope if and only if 

(9) conv(rows(M)) = aff(rows(M)) n R q + . 

Since we have 

conv(rows(M)) = <£>(conv(rows(V))) = <p(Q) 

and 

aff(rows(M)) nl| = ^(aif(rows(F))) DR q + = ip(aS(Q)) n R q + 

= p{{x G aff(Q) : <p(x) > O}) = <p(P n aff(Q)) = <p(P) 

(here we used that dim(P) = dim(Q)), condition ([9j is equivalent to ip(P) — <p(Q)- 
In turn, this is equivalent to P = Q since (p is injective. Thus, P — Q is equivalent 
to M being the slack matrix of a polytope. □ 

4. A Combinatorial Characterization of Slack Matrices 

Our second characterization of slack matrices of cones and polytopes relies on 
incidence structures. For a (nonnegative) matrix M, we denote by M; nc the 0/1- 
matrix with (-Mi nG )ij = 1 if and only if My = 0. The matrices M- mc arising from 
slack matrices M of a polyhedral cone if or of a polytope P are called the incidence 
matrices of K or P, respectively. 

We start by characterizing the slack matrices of polytopes, since the correspond- 
ing statement for cones can easily be deduced from the one for polytopes. The 
characterization is restricted to nonnegative matrices of rank at least two. It is 
easy to see that no matrix of rank one is a slack matrix of a nontrivial polytope. 
One may (or may not) want to consider a rank-zero matrix as a slack matrix of the 
polytope consisting of the zero- vector in R°. 

Theorem 22. A nonnegative matrix M with rank(M) > 2 is a slack matrix of 
some polytope if and only if M mc is an incidence matrix of some (rank(M) — 1)- 
dimensional polytope and 1 is contained in the column span of M . 

Proof. If M is a slack matrix of a polytope P, then 1 is contained in the column 
span of M (Theorem [6]) , and by Theorem [T4l dim(P) = rank(M) — 1. 
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In order to establish the non-trivial implication of the claim, let M £ K^. x<z be a 
nonnegative matrix with rank(M) = d + 1 > 2, 1 e M ■ 1' and M- lnc an incidence 
matrix of some d-dimensional polytope R. Denote by V C M.+ the set of rows 
of M and define the polytope P := conv(V / ) and the polyhedron Q := aff(V) n M. q + . 
Clearly, P C Q, and since 1 £ M ■ R q , dim(Q) = dim(P) = d. By Corollary in 
order to show that M is a slack matrix of a polytope, it suffices to prove P = Q. 

In order to establish Q C P, let us define 

14 = {t; G V : t?i = 0} and F t = conv(Vj) for 1 < i < q . 

The set 

i=l 

is contained in the relative boundary dQ of Q. Note that as an incidence matrix of 
some polytope of dimension at least one, Mj nc does not have an all-ones column. 
Since Q = conv(dQ) (note that Q is a pointed polyhedron of dimension d > 2, 
which is important here in case of Q being unbounded), if we show that F = dQ, 
then we will have that Q = conv(F) C P. 

Thus, our goal is to establish F = dQ. As mentioned above, we have F C dQ. 
It suffices to show that F is homotopy-equivalent to a (d — l)-dimensional spherqj, 
because then F cannot be properly contained in the (d — l)-dimensional connected 
(recall dim(Q) > 2) manifold dQ. This follows, e.g., from [3J Cor. 8.5] together 
with the fact that the (d — l)-st cohomology group of a (d — l)-dimensional sphere 
is non-trivial. 

To show that F is homotopy-equivalent to a (d— l)-dimensional sphere, observe 
that for every subset I C {1, , . . , g}, we have Hig/p ^ if and only if the submatrix 
of M; nc formed by the columns indexed by / has an all-ones row. Now let R be a 
polytope of which M- lnc is an incidence matrix. Let G\ , . . . , G q be the faces of R 
that correspond to the columns of M- mc . Then flig/Gj ^ holds if and only if the 
submatrix of M- mc formed by the columns indexed by / has an all-ones row. 

Therefore, the abstract simplicial complexes 

{IC{l,...,q}:f]Fi^0}, and {I C {1, . . . , q} : f| G, ^ 0} 

(known as the nerves of the polyhedral complexes induced by F\ , . . . , F q and by 
Gi, . . . , G q , respectively) are identical. Since all intersections Hie/ ft and f] ieI G{ 
are contractible (in fact, they are even convex), this simplicial complex is homotopy 
equivalent to both F and to the (d — l)-dimensional (polyhedral) sphere dR (see, 
e.g., |H Thm. 10.6]). ' □ 

Since polygons have a very simple combinatorial structure, Theorem readily 
yields a simple characterization of their slack-matrices. Here, a vertex-facet slack 
matrix of a polytope P is a slack matrix of P whose rows and columns are in 
one-to-one correspondence with the vertices and facets of P, respectively. 

Corollary 23. A matrix M G K" xn (n > 3 ) is a vertex-facet slack matrix of an 
n-gon if and only if its rows span an affine space of dimension exactly two and its 
rows and columns can be permuted such that the non-zero entries appear exactly at 
the positions (i,i) (for 1 < i < n), and (i,i — 1) (for 2 < i < n), and (l,n). 



Our proof of this is inspired by [7] . 
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Steinitz' theorem Q3| says that a graph G is the 1-skeleton of a three-dimensional 
polytope if and only if G is planar and three-connected. Using this, one can check 
in polynomial time whether a given 0/1-matrix is an incidence matrix of a three- 
dimensional polytope. For every fixed d > 4, however, it is NP-hard to decide 
whether a given / 1-matrix is an incidence matrix of a d-dimensional polytope |10) . 

In the following combinatorial characterization of slack matrices of cones we 
restrict our attention to matrices of rank at least two as for polytopes. Clearly, 
every nonnegative matrix of rank one is a slack matrix of the ray , and, we may 
consider a matrix of rank zero as a slack matrix of the trivial cone {0} in Mr. 

Theorem 24. A nonnegative matrix M with rank(M) > 2 is a slack matrix of 
a polyhedral cone if and only if M; nc is an incidence matrix of some rank(Af)- 
dimensional pointed polyhedral cone. 

Proof. If M is a slack matrix of some polyhedral cone then, by Lemma [TU1 M is 
a slack matrix (and hence M- lnc is an incidence matrix) of a pointed polyhedral 
cone K. By Lemma [T3l this cone has dimension rank(Af). 

In order to prove the reverse implication, we can assume by the results in Section 
2.3 that M does not have any zero-row. Since M is also nonnegative, there exists 
a positive diagonal matrix D such that DM contains 1 in its column span. 

Given a pointed cone K, we can slice K by an affine hyperplane L such that the 
slice is a polytope of dimension dim(_ftT) — 1 and the incidence structures of K and 
K fl L are identical. Thus, {DM), is an incidence matrix of some (rank(M) — 
l)-dimensional polytope. By Theorem [22j we have that DM is a slack matrix 
of a polytope. Hence, M is a slack matrix of the homogenization cone of this 
polytope. □ 

Note that dropping pointed from the formulation of Theorem [24] makes the state- 
ment false. Indeed, 



M 



1 2 

2 1 





With Mine = 







1 1 
1 1 



and rank(M) = 2 is not a slack matrix (since M does not satisfy the RCGC), 
but Mj nc is the incidence matrix of the non-pointed cone {(x\,X2) ■ x-i > 0} with 
V-representation (0,1), (0,1), (1,0), (—1,0) and "H-representation (0,1), (0,1). 
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